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ABSTRACT 

An intriguing hypothesis, that there exists a decaying X-ray emission component be- 
fore the GRB trigger, has been suggested in order to explain the shallow decay phase 
of the X-ray afterglow detected in many Swift GRBs. If this "prior emission" is from an 
external shock, one would expect a corresponding optical emission component during 
the GRB prompt emission phase. In this paper we apply the available prompt optical 
emission data (both detections and upper limits) to constrain such a scenario. We fit 
the shallow and normal decay segments of the XRT light curves in our sample with 
a TA-shifted single power law, and extrapolate the X-ray flux back to the time of the 
early optical observations. We then use the synchrotron spectrum predicted by the 
standard external shock model to extrapolate from the X-ray flux to the optical band, 
and obtain the possible range of the predicted optical flux. Finally, we compare the 
predictions with the observations. In the cases where later optical data are available, 
we also compare the shapes of the optical light curves to the predicted optical light 
curves from the external shock prior emission model. We find that for a good fraction 
of GRBs (4 out of 8; up to December 2006), the available data already impose severe 
constraints on the hypothesis. In particular, the expected optical flux from the prior 
external shock model is higher than what the data allow. We conclude that if the 
shallow-decay X-ray component were from a prior emission component, it would have 
to be of an internal origin with optical flux suppressed. 

Key words: gamma-rays: bursts, shock waves, radiation mechanisms: non-thermal 



1 INTRODUCTION 

Thanks to the vast collection of data on GRBs and their 
afterglows in the Swift era (Nousek et al. 2006; O'Brien 
et al. 2006; Willingale et al. 2007; Evans et al. 2009), we 
now have a detailed picture of the temporal behavior of the 
X-ray emission from these violent cosmic explosions. Using 
light curves obtained by the X-Ray Telescope (XRT) on- 
board the Swift satellite, five phases that are common to 
many X-ray afterglow light curves were identified: a steep 
decay phase, a shallow decay phase, a normal decay phase, 
a post jet break phase, and X-ray flares (Zhang et al. 2006). 
Not all bursts exhibit all five phases. The normal decay and 
post jet break phases were predicted pre-Swift and are well- 
explained by the standard external shock afterglow model 
(Meszaros & Rees 1997; Sari et al. 1998; Rhoads 1999; Sari et 
al. 1999). Swift observations have challenged the GRB com- 
munity to create physical models that account for the other 
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three phases. Of these phases, the shallow decay phase has 
definitely posed the biggest challenge. Generally, no spec- 
tral evolution is observed across the temporal break from 
the shallow to the normal decay phase (Liang et al. 2007). 
This suggests that the shallow-to-normal transition should 
be either hydrodynamical or geometrical in origin. Many 
such physical models of the shallow decay phase have been 
discussed in the literature, but none has been successful in 
interpreting all the observations (Zhang 2007 and references 
therein). 

An intriguing hypothesis was proposed by Yamazaki 
(2009). According to this hypothesis, there was already ac- 
tivity from the central engine before the GRB itself. The 
X-ray flux is already decaying before the GRB trigger. Due 
to a bad choice of reference time (the GRB trigger time), 
this power-law-decaying prior emission component (ejected 
at a time — Ta with respect to the GRB trigger time, which 
we define as T — 0) appears as an artificial broken power 
law with a plateau stretching out to T « Ta in the log-log 
diagram. In order to observe the single power law behavior 
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of the prior emission component, one needs to move the time 
zero point to — Ta, which we define as t — 0. According to 
Kobayashi & Zhang (2007), the time zero point of external 
shock emission should coincide with the beginning of the 
corresponding central engine activity. So for prior emission 
it should be at t = 0, while for the main GRB it should 
be at T = cQ. The observed X-ray afterglow is dominated 
by the prior emission component, which outshines the X- 
ray emission associated with the main GRB outflow during 
the shallow and normal decay phases. A systematic study 
of Swift/yLRT light curves from GRBs detected before mid- 
2009 (Liang et al. 2009) suggests that the prior emission 
model provides a unified interpretation of both the canon- 
ical X-ray light curves and the light curves that seem to 
follow a single power law decay (which are the ones with 
negligibly short Ta). In addition, Liang et al. (2009) also 
found that both types of light curves in their sample can be 
roughly explained by the external shock model. 

If a prior explosion indeed happened, this prior outflow 
must sweep up the ambient medium and form an external 
shock. The long-lasting prior X-ray emission that is invoked 
to interpret the X-ray plateaus is most likely from this ex- 
ternal shock. If this is the case, then the same external shock 
would give rise to an optical emission component that fol- 
lows the same plateau behavior. This can be tested with the 
available optical data, especially during the earliest obser- 
vational epoch. In general, optical emission observed during 
the prompt phase may be a mix of this prior emission com- 
ponent and several other emission components, including the 
optical counterpart of the prompt GRB emission, and the 
reverse shock and forward shock emission associated with 
the prompt outflow. As discussed in Ioka et al. (2006), the 
ejecta associated with the prompt emission would stream 
into the trail of the prior component blastwave with modi- 
fied medium density profile, leading to another pair of for- 
ward and reverse shocks (for the dynamics setting of the 
three shock system, see Zhang & Meszaros 2002). The rela- 
tive importance of emission contribution to the optical band 
depends on many parameters. The reverse shock component 
can dominate or be outshone by the forward shock compo- 
nent in the optical band. The detected prompt optical emis- 
sion flux is therefore allowed to be higher than the predicted 
flux. On the other hand, if the prompt optical flux or upper 
limit is already below the range of the predicted optical flux, 
then the external shock prior emission model would be ruled 
out. This is the motivation and strategy of this paper. 



2 SAMPLE SELECTION AND DATA 
REDUCTION 

In order to create a sample for this study, we need bursts 
with both X-ray plateaus (which can be fit within the prior 



1 Quimby et al. (2006) have explored the effects that changing 
the afterglow reference time has on the X-ray and optical light 
curves for the ambiguous case of GRB 050319, a burst in our 
sample. In our analysis of GRB 050319, we adopted their defined 
ttrO as T = 0, which marks the start of the prompt emission even 
though the actual trigger happened later due to an instrumental 
slewing effect. 



emission model) and early optical observations (either detec- 
tions or upper limits). We searched through the literature 
to find the GRBs which either clearly exhibit a shallow de- 
cay segment in their XRT light curve or have been fit with 
a shallow decay segment. Specifically, the candidate bursts 
were taken from Table 1 of Liang et al. (2007) and Table 2 
of Liang et al. (2009). We further narrowed down the sample 
to those GRBs with early optical observations, using Table 5 
from Yost et al. (2007a) and Table 2 from Yost et al. (2007b) 
for detections, and Table 3 from Yost et al. (2007b) for up- 
per limits. Thus, the effective cutoff date for our sample is 
December 2006, about two years into the Swift mission. Al- 
together, our sample consists of eight bursts (Table 1) with 
a total of 11 early optical detections and 17 early optical 
upper limits. In the interest of comparing the optical light 
curves of these GRBs to those predicted by the external 
shock prior emission model, we have included as much of 
the available optical data for these bursts as possible. The 
references for the optical light curve data are given in Ta- 
ble 1. Since most of the prompt optical observations from 
Yost et al. (2007a,b) are in the 7?-band, we include only 
the ij-band optical light curves. The only exception is GRB 
060927, for which the included prompt optical observations 
and optical light curve data are at a wavelength (8190 A) 
near the i-band (Ruiz-Velasco et al. 2007). 

Before we can compare the prompt optical observations 
and optical light curves from the literature to the predictions 
of the external shock prior emission model, the optical data 
must be corrected for Galactic and host galaxy extinction. 
We corrected the optical data for Galactic extinction using 
the empirical Milky Way (MW) extinction law, extinction 
coefficients, and Rv = 3.08 from Pei (1992), along with the 
values of As from Schlegel et al. (1998). In order to cor- 
rect the optical data for host galaxy extinction, we searched 
the literature on each burst in our sample for the extinction 
law and Av or Eb-v value which provided the best fit to 
its afterglow spectral energy distribution (SED). We then 
used the appropriate results from the literature to calculate 
the host galaxy extinction, once again relying on the empir- 
ical MW, Large Magellanic Cloud (LMC), and Small Mag- 
ellanic Cloud (SMC) extinction laws and parameters from 
Pei (1992). 

The host galaxy extinction information and references 
for each GRB are given in Table 1. For GRB 060927 and 
GRB 061222A, the cited authors were unable to distinguish 
between the three extinction laws. We chose the SMC ex- 
tinction law based on the conclusions of Schady et al. (2007, 
2010), that, in most cases, it provides an acceptable fit to the 
host galaxy extinction profile. It should also be noted that 
the host galaxy extinction of GRB 050401 is the subject of 
debate. We have decided to use the analysis of Watson et al. 
(2006), although higher values of Ay have been advocated 
by De Pasquale et al. (2006) and Kamble et al. (2009). 

The calculated values (in magnitudes) of the host 
galaxy extinction A\ emit and Galactic extinction A\ obs are 
given in Table 1. The error on A\ emit is due to propagat- 
ing the error on the host galaxy value of Av or Eb-v taken 
from the literature, and the error on Rv given by Pei (1992). 
Since the host galaxy A v values for GRB 051109A and GRB 
061222A are limits, the errors on their A\ cmit only take into 
account the error on Rv- We were unable to calculate the 
host galaxy extinction for GRB 050822 since there is no 
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Table 1. Optical light curve references, redshifts, host galaxy extinction information, values of host galaxy and Galactic extinction, and 
spectral fit results for our sample GRBs. 



GRB 


z 


Host galaxy 


Host 


Galactic 


Fitting 


Galactic 


Host 


P 


x 2 M 






extinction 


(mag) 




interval 


N H * 


N H * 










information 


(mag) 


(ks) 


(10 22 cm) 


(10 22 cm) 






050319" 


3.24 6 


A v = 0.21 ±0.08, MW C 


0.53 ±0.20 


0.028 


6.24-82.77 


0.0126 


0.34 


1.19 ±0.07 


0.83 (374) 


050401 d 


2.90 6 


A v = 0.62 ± 0.06, SMC e 


2.4 ±0.2 


0.167 


0.46-20.09 


0.044 


1.57 


1.00 ±0.07 


0.78 (365) 


050822 


1.2/ 


unknown 


? 




6.41-516.15 


0.014 


0.118 


0.90 ±0.04 


1.03 (421) 


051109A» 


2.35 6 


A v < 0.10, SMC C 


< 0.3189 ±0.0008 


0.487 


3.73-636.40 


0.16 


0.93 


1.29 ±0.03 


1.03 (560) 


060729 


0.54 b 


A v = 0.04 ±0.02, LMC C 


0.05 ±0.03 




0.70-2,221.24 


0.045 


0.099 


1.14 ±0.01 


1.16 (749) 


060927 h 


5.47 h 


E B -v = 0.07 ± 0.03 h , SMC 


1.1 ± 0.5 




0.11-10.86 


0.046 


0.860 


1.0 ±0.1 


0.50 (208) 


061121* 


1.31 6 


A v = 0.40 ±0.04, LMC C 


0.80 ±0.08 


0.116 


0.59-348.63 


0.04 


0.607 


1.09 ± 0.03 


0.88 (587) 


061222A 


2.088-? 


A v > 5.0?, SMC 


> 14.74 ±0.04 




0.23-723.42 


0.09 


3.82 


1.13 ± 0.03 


1.15 (692) 



*Taken from Evans et al. (2009) and held constant in spectral fit; "Optical light curve data from Greco et al. (2005a), Kiziloglu et al. 
(2005), Misra et al. (2005a), Sharapov et al. (2005), Wozniak et al. (2005), Quimby et al. (2006), Huang et al. (2007), and Kamble et al. 
(2007); b Liang ct al. (2007); c Schady et al. (2010); d Optical light curve data from D'Avanzo et al. (2005), Greco et al. (2005b), Kahharov 
et al. (2005), Misra et al. (2005b), Rykoff et al. (2005), De Pasquale et al. (2006), Watson et al. (2006), and Ghirlanda et al. (2007); 
e Watson et al. (2006); ^Butler (2007); ^Optical light curve data from Yost ct al. (2007a); ^Optical light curve data, redshift, and host 
galaxy extinction information from Ruiz-Velasco et al. (2007); 'Optical light curve data from Efimov et al. (2006), Halpern et al. (2006), 
Halpern & Armstrong (2006a,b), Melandri et al. (2006), Page et al. (2007), and Uehara et al. (2011); ^Perley et al. (2009) 



information about its host galaxy extinction profile in the 
literature. We did not calculate Galactic extinction for GRB 
050822, GRB 060729, GRB 060927, and GRB 061222A since 
all of the included optical flux densities for these bursts have 
already been corrected for Galactic extinction by either Yost 
et al. (2007a,b) or Ruiz-Velasco et al. (2007). The errors on 
A\ ob3 were negligible. 

The Swift/XRT data for the GRBs in our sample are 
taken from the Swift data archive^. We developed an IDL 
package which invokes the standard HEASoft tools (e.g., Xs- 
elect, Ximage, Xspec, Swift tools) to automatically process 
the XRT data for any given burst. The details of the proce- 
dure of our package are described in Zhang et al. (2007b). 
Our final data product^] include the XRT light curves (in 
physical units), extracted over the energy range 0.3 — 10 keV. 
These data products have been applied to study statistical 
properties of X-ray afterglows in a series of papers (Zhang 
et al. 2007b; Liang et al. 2007, 2008, 2009). 

For the bursts in our sample, the spectral index /3 for the 
XRT band is obtained by fitting the time-integrated spec- 
trum (Evans et al. 2009) of the shallow and normal decay 
phases with a simple power law model in Xspec. The spec- 
tral fit also includes two-component neutral hydrogen col- 
umn density Nn.host and Nn.Gai, with corresponding val- 
ues taken from the late-time spectral data in the Swift/XRT 
GRB spectrum repository (Evans et al. 2009). Their values 
were held constant in our spectral fit, which are collected 
in Table 1 along with the GRB redshift z and the spectral 
fitting results. 

All the errors quoted in this paper correspond to the 
la errors on the quantities, unless otherwise indicated. The 
convention F u cc t' - a v-P is adopted. The zero points and 
effective frequencies for the optical bands from Bessell (1979) 
are adopted. 



2 ftp : //legacy . gsf c .nasa. gov/swif t/data/ 

3 Available from http://grb.physics.unlv.edu/swift 

4 http : / /www . swift . ac . uk/xrt_spec tra/ | 



3 MODEL 

In this paper we focus on the XRT data in the shallow and 
normal decay phases (dark black points in Fig. 3). We do not 
consider the post jet break phase as the jet break happens 
much later. Using the prior emission model and Eq. [2] from 
Liang et al. (2009) as our prescription, we fifl the shallow 
and normal decay segments of each XRT light curve to the 
function 

allowing the normalization constant Fo, time shift constant 
Ta, and temporal decay index a to vary in the fit. Here 
T is the time measured since the GRB trigger. The fitting 
intervals are loosely based on the time ranges given in Ta- 
ble 1 of Liang et al. (2007). Once the parameters have been 
determined for a particular burst, Eq. [TJ can be used to 
extrapolate the X-ray flux due to the prior emission compo- 
nent from the shallow decay segment to earlier times (i.e., to 
the time of the early optical observations of interest). The 
results of the XRT light curve fits are given in Table 2. 

The external shock model predicts a broken power law 
synchrotron spectrum (Meszaros & Rees 1997; Sari et al. 
1998). For a constant density medium, a fast cooling spec- 
trum (y c < Vm) is expected at early times, which transitions 
to a slow cooling spectrum (v m < v c ) at later times (Sari et 
al. 1998). In our analysis, since T A ~ 10 3 - 10 4 s (typical 
observed shallow-to- normal break time), the prior emission 
external shock should have already entered the slow cooling 
regime by the time of the GRB trigger. Thus, we apply slow 
cooling throughout the modeling. 

We first determine the spectral regime in which the 
XRT band is located by using the X-ray a and /3 data and 
comparing the data with the closure relations of the exter- 
nal shock models (Zhang & Meszaros 2004). Our approach 
is similar to that used by Racusin et al. (2009). We use the 

5 Using the routine MPFITFUN.PRO (Markwardt 2009): 
http : / /purl . com/net/mpf it 
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Table 2. Results of temporal fits and spectral regime/model determinations. 



GRB Fitting F T A a \ 2 { v ) Regime Model ot pred 

interval (10 -11 erg (ks) 

(ks) s _1 cm 2 ) 



050319 


6.11-84.79 


2.1 ± 0.3 


17 ±10 


1.5 ±0.4 


0.51 (27) 


vx 


> v c 


ISM, p > 2 


1.29 ±0.11 


050401 


0.14-20.09 


87 ±4 


0.53 ± 0.05 


1.13 ±0.03 


2.21 (132) 


vx 


> V c 


ISM, p > 2* 


1.00 ±0.11 


050822 


6.41-523.32 


1.5 ± 0.3 


12 ±5 


1.1 ±0.1 


1.24 (72) 


v m 


<vx <v c 


ISM, p > 2* 


1.35 ±0.06 


051109A 


3.73-639.16 


15 ±2 


3.3 ±0.7 


1.33 ±0.04 


1.17 (94) 


vx 


> V c 


ISM, p > 2* 


1.44 ±0.05 


060729 


0.70-2,221.24 


2.93 ±0.05 


120 ±7 


1.96 ±0.05 


1.69 (813) 


v m 


< v x < v c 


ISM, p > 2* 


1.71 ±0.02 


060927 


0.11-10.86 


4.4 ±0.5 


2.7 ±1.8 


2.5 ± 1.0 


1.19 (29) 


Vm 


<Vx<V c 


Wind, p > 2 


2.0 ±0.2 


061121 


0.30-353.10 


35.8 ±0.9 


2.6 ±0.1 


1.35 ±0.02 


1.18 (412) 


vx 


> V c 


ISM, p > 2* 


1.14 ±0.05 


061222A 


0.23-724.62 


10.2 ±0.2 


8.5 ±0.4 


1.56 ±0.02 


1.72 (317) 


Vm 


<v x <v c 


ISM, p > 2* 


1.70 ±0.05 



*a and ct pre d not consistent within error for this regime and model 



value of (3 from our spectral fit to evaluate the slow cool- 
ing closure relations, with each model/regime combination 
yielding a predicted value of the temporal decay index a pre d- 
We then compare the values of a pr ed to the value of a ob- 
tained from our temporal fit, looking for the model/regime 
combination for which a — a pre d = (within lcr error). Of 
course, the decision typically was not clear cut. For six out 
of the eight GRBs in our sample, a was not consistent with 
ctpred from any of the model/regime combinations within 
error. In these cases, we chose the model and regime with 
which the data was most consistent (minimizing \a — a pre d\)- 

Ultimately, the XRT band is found to be either in the 
v m < vx < v c or the vx > v c spectral regimfl The results 
of our spectral regime and model determinations are given 
in Table 2. One ambiguous case was GRB 050822. Its XRT 
data were nearly equally consistent with three model/regime 
combinations: vx > v c for ISM and Wind [p < 2), and v m < 
vx < v c for ISM (p > 2). We chose the last model/regime 
combination for GRB 050822 because the resulting predicted 
range for the optical flux actually encompasses the predicted 
ranges corresponding to the other two model/regime choices. 

Next, we apply the XRT light curve fit results to ex- 
trapolate the prior emission X-ray flux back to the time of 
the optical observations (T ob s), and derive the mean X-ray 
spectral density as 



Fx{T ob 



vx 



T bs + Ta 



(2) 



where vx « 4.8 x 10 17 Hz (2 keVfl is applied. 

We then apply the slow cooling model (Fig. lb of Sari et 
al. 1998) to predict the range of optical flux density. In the 
case that vx > v c (Fig. la), the maximum possible optical 
flux density occurs when v opt > v c , so that the optical band 
lies in the same spectral regime and one can do a simple 
extrapolation of X-ray flux density based on the power law 
spectral fits (top dashed line): 



Fopt,max(F bs^) — FxiTobs') 



( Vopt \ 

V vx J 



(3) 



6 In some cases, a burst is consistent with both ISM and Wind 
models within the same spectral regime. For the analysis below, 
the medium density profile does not enter the problem. So we 
care about the spectral regime only. 

7 Approximate log average of the XRT energy range 



If v op t is below v c , then most likely v m < v op t < v c . The 
minimum possible optical flux density in this regime occurs 
when vx is just above v c , so that 



F 



opt,min 



(T bs) — Fx (T bs) 



( Vopt 
V V X 



-(/3-1/2) 



(4) 



This is marked as the second dashed line in Fig. la. In the 
most extreme cases, one may have both v m and v c between 
the optical and X-ray bands. Therefore the lowest flux den- 
sity allowed corresponds to v m — v c < vx, so that 



F pt,min(T bs) — Fx (T bs) 



( Vopt \ 

V vx J 



1/3 



(5) 



This corresponds to the third dashed line in Fig. la. 

In the case that v m < vx < v c (Fig. lb), the maximum 
possible optical flux density occurs when vx is just below v c 
and Vopt is in the same regime as vx (Eq. [3]). The minimum 
possible optical flux occurs when vx is just above v m and 
Vopt < v m < v x (Eq. [5]). 

In both Fig. la and Fig. lb, the thicker part of the 
range corresponds to the maximum prediction, for which 
there are no spectral breaks between the X-ray and optical 
bands. The thinner part of the range, which includes the 
minimum prediction, represents the case in which there is 
one spectral break between the bands. The dash-dotted line 
in Fig. la corresponds to the extreme case in which there 
are two spectral breaks between the X-ray and optical bands. 
In all cases, we have introduced an error on the prediction, 
which includes the uncertainties in ft, Fq, Ta, and a, with 
the prior emission model and external shock synchrotron 
model applied and errors properly propagated. 

The final predicted range for the optical flux density can 
then be compared to the optical observations with proper ex- 
tinction correction. If all of the extinction-corrected prompt 
optical observations fall within their predicted ranges, then 
the data would show no inconsistency with the hypothesis of 
an external-shock-origin prior emission. If some or all of the 
observations lie above their predicted ranges, this would not 
count as evidence against this hypothesis, because the ob- 
served optical flux can have contributions from the internal 
dissipation region and the reverse/forward shock associated 
with the prompt emission, besides the contribution from the 
prior emission component. However, an observation that lies 
below its predicted range would put a severe constraint on 
the external shock prior emission model. 
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Optical 



XRT 



logv 



> 

v 



P~-(p-l)/2 



p~ 1/3 



Optical 



XRT 



logv 



Figure 1. Cartoon spectra based on the results of Sari et al. 
(1998). For each possible spectral regime for i/j, there are min- 
imum and maximum possibilities for the corresponding optical 
flux density (section 3). 



0.01 



1000 
Time (s) 



10' 



L0 5 



Figure 2. Example optical (solid) and X-ray (dashed) light 
curves derived from the external shock prior emission model 
with typical parameter values: t m = 2500 s (p = 2.3, e e = 0.3, 
e s = 0.01, E K = 10 52 erg, n = 0.5, z = 1.0) and T A = s 
(black), 500 s (red), 1250 s (green), 2500 s (blue), 10000 s (pink). 
These light curves were generated by our GRB afterglow code 
which uses the results of Sari et al. (1998) and Zhang et al. (2007a) 
as a prescription. 



The interpretation becomes more complicated in the 
cases where later optical observations are available. In these 
cases, one may examine the overall shape of the optical light 
curve to determine whether or not it is consistent with the 
external shock prior emission model. If the prior emission 
time is chosen as the zero point, the X-ray lightcurve and 
optical lightcurve typically have different shapes. In partic- 
ular, if v opt < v m early on, the optical lightcurve would 
display an initial rise with F opt oc t 1 ' 2 (where t is time mea- 
sured from the prior emission zero point) followed by a decay 
(Sari et al. 1998). With T A shift, the shapes of both X-ray 
and optical lightcurves are distorted. In particular, the op- 
tical lightcurve as measured since the GRB trigger may be 
well fit by a smooth broken power law 



F = F 



T + T A 

T m + Xa 



+ 



T + T A 

T m + Ta 



(6) 



where T m is the time referenced from the GRB trigger when 
v m crosses v opt and w is the smoothing factor (Liang et al. 
2007). Figure 2 shows the comparison between the X-ray 
and optical lightcurves for various values of Ta- It is found 
that the difference between the shapes of the two lightcurves 
reduces as Ta increases, and essentially disappears as Ta 
exceeds t m (the time referenced from prior emission when 
v m crosses Vopt). 

If the optical lightcurve does not follow the predicted 
shape of the external shock prior emission model, one would 
conclude that it is dominated by another component. In 
that case, the allowed optical emission level from the ex- 
ternal shock prior emission component should be at most 
the highest predicted prior emission lightcurve allowed by 
the data. Additionally, the predicted prior emission optical 
light curve should lie within the predicted range delineated 
by the bars in Figure 1. 
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10-* 10* 10 
Time since trigger (s) 



10' 10 4 10' 
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Figure 3. The grey and black asterisks are the XRT light curve data. The red line is the fit of Eq. [I] to the black XRT light curve data 
points. The green bars represent the predicted ranges for the optical flux density at T >, s . The thicker part of the green bar corresponds 
to the maximum prediction for which there are no spectral breaks between the XRT and optical bands. The thinner part of the green 
bar, which includes the minimum prediction, is obtained by assuming there is a spectral break between the bands. For GRB 061121, the 
dash-dotted portion of the green bar represents the extreme minimum prediction obtained by assuming there are two spectral breaks 
between the XRT and optical bands. The dark blue circles and triangles represent the optical detections and upper limits, respectively. 
The dashed black lines are the predicted optical light curves from the external shock prior emission model, as constrained by the green 
ranges and the optical data. For GRB 061121, the dashed light blue line is the fit of Eq. [l] to the optical light curve data in which 
is fixed to the value from the XRT light curve fit. 
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4 RESULTS AND DISCUSSION 

The results of the spectral fits, burst redshifts, extinction in- 
formation, and the inclusive list of the references for the opti- 
cal light curve data are given in Table 1. The XRT light curve 
fits and spectral regime determinations are given in Table 
2. The XRT light curve data, associated fits, extinction- 
corrected optical observations, and predicted ranges are 
shown in Figure 3. The green bars in the plots have the 
same meanings as the bars shown in the cartoon spectra in 
Figure 1 and explained in section 3. The dashed black lines 
are the maximum predicted optical light curves from the 
external shock prior emission model allowed by the data. 
If the optical data constrain the dashed black line to be 
below the green range, for self-consistency we instead plot 
a dashed line corresponding to the minimum optical light 
curve allowed by the external shock prior emission model, 
even though it clearly contradicts the data. The dashed light 
blue line in the plot for GRB 061121 is a fit to its optical 
light curve data (discussed in section 4.2). Looking at Fig- 
ure 3, one can see that the results for each GRB are a bit 
different. In general, we can assign the bursts to one of two 
groups: those that are consistent (or at least, not inconsis- 
tent) with the external shock prior emission hypothesis, and 
those that are inconsistent with the external shock prior 
emission hypothesis. 



4.1 "Consistent" Bursts: GRB 050822, GRB 
060729, GRB 060927, and GRB 061222A 

The early optical upper limits of GRB 050822 lie within the 
thick part of the green range, corresponding to the maximum 
prediction in which v m < v opt < vx < v . However, since 
these data are just upper limits and not actual detections, 
the actual level of optical flux from this dark burst could lie 
anywhere in the green range or even below the green range. 
In fact, since we were unable to correct these upper limits 
for host galaxy extinction, the limits (and thus the actual 
level of optical flux) could even be above the green range. 
Regardless, the facts that these are only upper limits, and 
that even without correction for host galaxy extinction they 
lie within the green range, lead us to the conclusion that 
GRB 050822 is not inconsistent with the external shock prior 
emission hypothesis. 

GRB 060729 has been studied (e.g., by Grupe et al. 
2007, Curran et al. 2009) in many optical bands. However, 
there is no additional 7?-band data in the literature. We 
find the early i?-band optical upper limit and detections of 
GRB 060729, which have been corrected for both Galactic 
and host galaxy extinction, lie within the thinner part of 
the green range. This is consistent with the case in which 
there is one spectral break between the X-ray and optical 
bands: v opt < v m < fx < v c . Due to the lack of additional 
_R-band data, we cannot make any further comparisons of 
this burst to predictions of the external shock prior emission 
hypothesis, and conclude that the available 7?-band data for 
GRB 060729 are not inconsistent with this hypothesis. 

For GRB 060927, the early optical data is at a wave- 
length near the i-band, and fortunately there is additional 
optical data in this band in the literature. The prompt «- 
band optical flux density falls within the thicker part of the 
green predicted range, corresponding to the maximum pre- 



diction in which v m < v op t < vx < ix c . It is clear that the 
overall shape of the optical light curve is not consistent with 
the smooth, TA-shifted single power law predicted by the 
external shock prior emission model. This suggests that, if 
there is a prior emission external shock optical component, 
it is dominated by other optical components. The dashed 
black line shows the maximum prior emission optical light 
curve allowed by the data. 

The results for GRB 061222A are similar to those for 
GRB 050822, except in this case, the early optical upper 
limits are nearly four orders of magnitude above the green 
range. This is explained by the large amount of host galaxy 
extinction (nearly 15 magnitudes) for this burst. Because 
GRB 061222A is a dark burst, we don't know whether the 
actual level of optical flux is consistent with the prediction 
of the external shock prior emission model. Thus, based on 
the limited data available, we conclude that this burst is not 
inconsistent with our hypothesis. 

Overall, the optical data of GRB 050822, GRB 060729, 
GRB 060927, and GRB 061222A are not inconsistent with 
the hypothesis that the prior emission, if it exists, arises 
from an external shock. However, due to the lack of optical 
light curve data for comparison (in the cases of GRB 050822, 
GRB 060729, and GRB 061222A), or excess optical flux (in 
the case of GRB 060927), these bursts do not make a strong 
case for the hypothesis. In other words, the results for these 
GRBs do not constrain or rule out the external shock prior 
emission hypothesis. 

4.2 Inconsistent Bursts: GRB 050319, GRB 
050401, GRB 051109A, and GRB 061121 

The early optical detection of GRB 050319 falls within the 
thinner part of the green range. This is consistent with the 
prediction in which there is one spectral break between the 
X-ray and optical bands: v m < v op t < v c < vx- How- 
ever, the optical light curve data clearly do not exhibit the 
smooth, TA-shifted single power law behavior predicted by 
the external shock prior emission model. This suggests that, 
even if the prior emission does originate from an external 
shock, the optical prior emission component is dominated 
by other optical components. Thus, the maximum predicted 
prior emission optical light curve should lie below the ac- 
tual optical light curve data. However, plotting the dashed 
line below all of the optical data would place the predicted 
prior emission light curve below the green range, which is 
not self-consistent. Thus, the dashed black line is ultimately 
constrained by the green range in this case. We must con- 
clude that GRB 050319 is not consistent with the external 
shock prior emission hypothesis. 

The situation is similar for GRB 050401 and GRB 
051109A. Although the early optical detections of these 
bursts fall within the thinner part of the predicted green 
range corresponding to v m < v op t < v c < vx, the shapes 
of the optical light curves are not consistent with the shape 
predicted by the external shock prior emission model. Once 
again, the predicted prior emission optical light curves are 
ultimately constrained by the bottom of the green range, 
and we must conclude that GRB 050401 and GRB 051109A 
are not consistent with the external shock prior emission 
hypothesis. 

GRB 061121 is the only burst in our sample for which 
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the early optical data fall below the green range. Conse- 
quently, the only way to explain this burst within the exter- 
nal shock prior emission model is to introduce the ad hoc 
assumption that there are two spectral breaks between the 
X-ray and optical bands {v opt < v m < v c < Vx), as indi- 
cated by the dash-dotted extension of the green range. Even 
after making this ad hoc assumption, there are still issues 
with GRB 061121. At first glance the shape of the optical 
light curve appears promisingly similar to that of the X-ray 
light curve. This could be consistent with Eq. [6] and Figure 
2 (given the right set of parameter values) , considering that 
Ta ~ 2500 s for this burst. However, when we actually fit 
the optical light curve (minus the second data point, which 
we dismissed as an optical flare) with Eq. [TJ, the results 
are not consistent with those of the X-ray light curve fit: 
T A = 1100 ± 300 s,a = 1.03 ± 0.02, x 2 0) = 0.83 (35). We 
then tried fitting the optical data, holding the time-shift 
parameter fixed at the value from the X-ray light curve 
fit (Ta = 2600 s), but the results still were inconsistent 
with those of the X-ray light curve fit: a — 1.09 ± 0.01, 
X 2 {v) = 1.03 (36) (this fit is shown as the dashed light 
blue line in Figure 3). Finally, we tried fitting the optical 
data with both the time-shift parameter and the decay in- 
dex fixed at the X-ray values (Ta = 2600 s, a = 1.35), but 
this yielded an unacceptable x 2 (y) = 11.2 (37). Due to the 
contrived parameter space required to make this burst con- 
sistent with the predictions, as well as the decay index of the 
optical light curve, we conclude that GRB 061121 is incon- 
sistent with the external shock prior emission hypothesis. 

Taken together, the optical data of GRB 050319, GRB 
050401, GRB 051109A, and GRB 061121 (which make up 
half of our sample), severely constrain the external shock 
prior emission hypothesis. 



5 CONCLUSION 

In order to test whether or not the observations are consis- 
tent with the intriguing hypothesis that the observed X-ray 
afterglow plateaus are due to an external shock prior emis- 
sion component, we have compared the prompt (and late, 
if available) optical emission data with the predicted ranges 
based on the model. We find that in 4 out of 8 GRBs, the 
available optical data already impose severe constraints on 
the hypothesis. In particular, the expected optical flux from 
the prior external shock model is higher than what the data 
allow. This suggests that the X-ray emission, if it has a prior 
emission origin, cannot be synchrotron emission from the ex- 
ternal shock. If indeed the shallow-decay X-ray component 
were from a prior emission component, it would have to be of 
an internal origin with optical flux suppressed. Future obser- 
vations before GRB triggers would definitely test the prior 
emission hypothesis. This might be realized in the SVOM 
era (see e.g. Paul et al. 2011 for a discussion of the SVOM 
mission and Zhang 2011 for a discussion of the prospects of 
detecting prior emission), or perhaps as soon as tomorrow, 
via serendipitous observations by current observatories. 
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